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Preface 


This  investigation  was  supported  by  the  Aquatic  Plant  Control  Research 
Program  (APCRP),  sponsored  by  the  Office,  Chief  of  Engineers  (OCE) ,  and  was 
managed  by  tha  US  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Hiss.  The  OCE  Technical  Monitor  was  Mr.  E.  Carl  Brown. 

This  is  the  final  report  for  Contract  No.  DACW39-81-C-0036,  "A  Mathemati¬ 
cal  Model  of  Submersed  Aquatic  Plants,"  prepared  by  Rensselaer  Polytechnic 
Institute  (RPI) ,  Troy,  N.  Y.  Authors  of  this  report  were  Drs.  Carol  Desormeau 
Collins,  Richard  A.  Park,  and  Charles  W.  Boylen,  RPI.  The  model  was  concep¬ 
tualized  and  developed  for  incorporation  into  the  US  Army  Corps  of  Engineers' 
reservoir  water  quality  model,  CE-QUAL-R1,  which  was  developed  during  the  con¬ 
duct  of  the  Environmental  and  Water  Quality  Operational  Studies  (EWQOS) . 
CE-QUAL-R1  is  a  numerical,  one-dimensional  model  that  describes  the  vertical 
distribution  of  thermal  energy  and  biological  and  chemical  materials  in  a  res¬ 
ervoir  through  time.  The  mathematical  structure  of  the  model  is  based  on  hori¬ 
zontal  layers;  temperature  and  materials  concentration  gradients  are  computed 
only  in  the  vertical  direction. 

The  original  contract  called  for  the  development  of  algorithms  and  the 
programming  of  those  algorithms  for  Inclusion  in  CE-QUAL-R1.  However,  in  sub¬ 
sequent  discussions  with  the  contract  officer  at  the  time,  Mr.  Joseph  Norton, 
Environmental  Research  and  Simulation  Division  (ERSD) ,  and  with  other  staff  of 
the  WES,  Environmental  Laboratory  (EL),  including  Drs.  Joseph  H,  Wlosinski  and 
Allan  S.  Lessem,  it  was  agreed  that  the  programming  should  be  done  by  the 
Environmental  Laboratory  staff  most  familiar  with  CE-QUAL-R1.  The  draft 
report  was  reviewed  by  Drs.  Wlosinski  and  Lessem  and  Messrs.  Mark  S.  Dortch 
and  Jack  B.  Walde. 

Manager  of  the  APCRP  was  Mr.  J.  Lewis  Decell.  General  supervision  was 
provided  by  Mr.  Donald  L.  Robey,  Chief,  ERSD.  Chief  of  the  EL  during  the  con¬ 
duct  of  this  investigation  was  Dr.  John  Harrison. 

Commanders  and  Directors  of  WES  during  the  study  and  preparation  of  the 
report  were  COL  Tilford  C.  Creel,  CE,  and  COL  Robert  C.  Lee,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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This  report  should  be  cited  as  follows: 

Collins,  C.  D.,  Park,  R.  A.,  and  Boylen,  C.  W.  1985.  "A 
Mathematical  Model  of  Submersed  Aquatic  Plants,"  Miscellaneous 
Paper  A-85-2,  prepared  by  Rensselaer  Polytechnic  Institute, 
Troy,  N.  Y.,  for  the  US  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  Miss. 
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A  MATHEMATICAL  MODEL  OF  SUBMERSED  AQUATIC  PLANTS 

Introduction 


Background 

1.  Submersed  aquatic  plants  or  macrophytes  often  contribute  signifi¬ 
cantly  to  the  productivity  of  lakes  and  reservoirs.  Macrophytes  can  become  so 
abundant  that  they  become  a  nuisance  to  recreational  and  navigational  activi¬ 
ties.  Their  growth  and  decomposition  also  influence  other  biotic  and  abiotic 
components  of  the  ecosystem.  The  littoral  community  of  many  eutrophic  systems 
is  often  dominated  by  a  single  species  of  macrophyte.  Under  less  eutrophic 
conditions,  'several  species  may  coexist.  The  growth  of  aquatic  plants  Is  con¬ 
trolled  by  many  factors,  including  (a)  growth  properties  of  the  plant; 

(b)  physical  factors  such  as  temperature,  irradiance  levels,  and  changes  in 
water  elevation;  and  (c)  physiological  characteristics  of  the  plant  such  as 
nutrient  requirements,  photoadaptation,  and  sediment  preference. 

2.  The  importance  of  macrophytes  to  the  aquatic  ecosystem  necessitated 
the  development  and  incorporation  of  a  macrophyte  submodel  in  the  US  Army 
Corps  of  Engineers'  one-dimensional  reservoir  water  quality  model,  CE-QUAL-R1 
(Environmental  Laboratory  1982) ,  which  was  developed  during  the  conduct  of  the 
Environmental  and  Water  Quality  Operational  Studies  (EWQOS) .  This  report 
describes  the  development  and  formulation  of  this  macrophyte  submodel  for 
inclusion  in  CE-QUAL-R1.  The  model  simulates  growth  and  decomposition  of 
macrophytes.  The  influence  of  the  plants  on  other  compartments  in  CE-QUAL-R1 
is  also  included  in  the  model. 

3.  To  make  the  proposed  submodel  complementary  with  CE-QUAL-R1,  the 
following  recommendations  are  made  regarding  the  computation  and  layering 
scheme  of  CE-QUAL-R1.  Macrophytes  should  be  regarded  as  occupying  the  bottom 
surface  of  each  layer  in  the  reservoir  within  the  euphotic  zone.  As  such, 
they  are  not  subject  to  advection  or  diffusion  and  are  not  transported  in  in¬ 
flowing  or  outflowing  waters.  The  macrophyte  compartment  should  have  units  of 
grams  per  layer.  As  the  layers  are  resized  in  CE-QUAL-R1,  dependent  on  the 
balance  of  inflowing  and  outflowing  waters,  the  macrophyte  biomass  should  be 
reapportioned  to  reflect  the  appropriate  densities  for  those  layers.  If  the 
surface  elevation  drops,  macrophytes  in  the  dewatered  zone  should  no  longer  be 
included  in  the  computation.  If  the  water  surface  elevation  increases  and 


inundates  new  areas,  the  macrophyte  density  in  the  new  area  should  be  given  a 
small  "seed"  value  to  represent  colonization. 

4.  Irradiance  reaching  a  particular  model  layer  determines  the  plants' 
growth  response.  Changes  in  water  level  can  affect  irradiance  at  a  particular 
level.  Drawdown  may  suddenly  expose  submersed  plants  to  higher  irradiances  as 
the  depth  of  water  through  which  light  is  transmitted  decreases.  Conversely, 
an  increase  in  reservoir  pool  elevation  may  result  in  greater  light  attenua¬ 
tion.  Light  attenuation  for  a  particular  layer  in  CE-QUAL-R1  is  dependent 
upon  the  extinction  coefficient  of  water  and  on  shading  by  suspended  solids, 
detritus,  zooplankton,  and  phytoplankton.  It  is  recommended  that  self-shading 
for  macrophytes  also  be  included  in  the  model. 

5.  The  following  processes  are  recommended  for  inclusion  in  the  macro¬ 
phyte  model:  gross  production,  dark  respiration,  photorespiration,  nonpreda- 
tory  mortality,  and  grazing.  Control  measures  affecting  macrophytes,  such  as 
mechanical  harvesting  and  herbicidal  treatment,  should  also  be  included  in  the 
model  as  described  in  this  report.  Decomposition  processes  already  modeled  in 
CE-QUAL-R1  would  be  affected  by  macrophyte  contributions  to  existing  detritus 
and  sediment  compartments.  A  flow  diagram  of  the  interactions  of  the  new 
macrophyte  compartment  with  other  model  compartments  summarizes  the  proposed 
changes  to  CE-QUAL-R1  (Figure  1). 


Figure  1.  Compartment  diagram  of  macrophyte  model  recommended  for  CE-QUAL-R1 
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Report  composition 

6.  In  the  following  section  the  specific  physiological  processes  recom¬ 
mended  for  inclusion  in  a  new  macrophyte  subroutine  are  formulated  for  incor¬ 
poration  into  CE-QUAL-R1.  Next,  a  geometric  scheme  for  apportioning  macrophyte 
biomass  among  model  layers  is  discussed.  The  next  major  section  contains 
recommendations  for  the  simulation  of  macrophyte  control  measures  (mechanical 
harvesting,  herbicidal  treatment).  The  next  section  discusses  the  validation 
of  select  process  formulations  based  upon  published  data  on  two  macrophyte 
species,  Myriophyllum  spicatum  and  Hydrilla  vertioillata .  The  final  section 
summarizes  the  major  recommendations  contained  in  this  report.  Two  appendices 
are  also  included.  Appendix  A  presents  equations  included  in  a  stand-alone 
version  of  the  macrophyte  submodel  used  in  the  process  validation  studies, 
while  Appendix  B  lists  representative  values  for  parameters  included  in  the 
proposed  macrophyte  submodel  based  on  published  research  on  M.  spicatum  and 
H.  vertioillata.  The  material  contained  in  this  report  will  be  included  in  a 
final,  revised  edition  of  the  CE-QUAL-RL  User's  Manual  (Environmental  Labora¬ 
tory  1982)  scheduled  for  publication  in  1985. 


Recommended  Physiologic  Processes 


7.  The  differential  equation  for  the  macrophyte  state  variable  expresses 
conservation  of  mass  in  each  horizontal  model  layer.  The  solution  provides 
material  concentrations  as  functions  of  time  and  depth.  The  equation  is  mathe¬ 
matically  expressed  as  follows: 


rate  of 
change 
of  mass 

g  day  1 


macrophyte 

biomass 


gross 

*  production 
rate 


dark 

respiration 

rate 


photorespiration 

rate 


nonpredatory 
grazing  -  mortality 
rate  rate 


mechanical 

or 

chemical 

harvesting 

rate 


(1) 


Each  of  the  individual  terms  in  this  equation  is  discussed  in  the  subsections 
which  follow.  The  style  of  presentation  follows  that  contained  in  the 


CE-QUAL-R1  User's  Manual  (Environmental  Laboratory  1982)  which  should  be  con¬ 
sulted  for  further  details.  The  overall  structure  of  CE-QUAL-R1  will  not  be 
presented  here.  Only  those  macrophyte  process  terms  specifically  included  in 
the  proposed  new  macrophyte  submodel  will  be  documented  plus  their  interac¬ 
tions  with  other  compartments  in  CE-QUAL-R1. 


Macrophyte  processes 


8.  Gross  production. 


The  daily  photosynthetic  or  gross  production  rate 


is  a  function  of  temperature,  light  intensity,  and  nutrient  concentration: 


PLTGRO  -  PLTMAX  *  RMULTl(T)  *  RMULT2(T)  *  MIN(XLIMN,XLIMP,XLIMC)  *  XLIML  (2) 


where 

PLTGRO  =  photosynthetic  rate,  day 

PLTMAX  *  user-specified  maximum  photosynthetic  rate,  day  * 
RMULT1,2(T)  =  temperature  limitation  functions,  unitless 
XLIMN  «  limitation  function  for  nitrogen,  unitless 
XL IMP  «  limitation  function  for  phosphorus,  unitless 
XLIMC  =  limitation  function  for  carbon,  unitless 
XLIML  -  limitation  function  for  light  intensity,  unitless 
9.  Temperature  limitation  is  calculated  using  the  equations  developed 
by  Thornton  and  Lessem  (1978): 


RMULTl(T)  = 


RMULT2(T)  - 


T  <  Ti 
T  >  Tl 

(3) 

T  <  T 

4 

T  £  T, 

4 


where 


K2(l  -  Kl) 


1 


X 


i 


X 


2 


In 


In 


Kjd  -  K2) 
K3(i  -  K4) 

V1  -  ty 


As  is  the  case  in  the  parent  model  CE-QUAL-R1,  and  represent  the 

user-specified  lower  and  upper  lethal  temperatures  for  the  processes  in 
question,  while  T2  and  (also  user  specified)  define  the  range  of 

optimum  temperatures  over  which  the  process  occurs  at  near  the  maximum  rate 
(Environmental  Laboratory  1982) .  The  term  T  represents  the  computed 
temperature  of  a  specific  layer  in  the  model  CE-QUAL-R1.  The  corresponding 
user-specified  K  values  define  the  relative  rates  (i.e.,  on  a  0  to  1  basis) 
at  which  the  process  occurs  at  each  of  these  temperatures. 

10.  Nutrient  limitation  is  dependent  upon  the  concentrations  of  nitrogen 
and  phosphorus  in  the  water  column  and  sediment  and  on  the  carbon  concentration 
in  the  water  column.  The  nutrient  determined  to  be  limiting  based  upon  the 
following  Monod  equation  is  used  in  the  photosynthesis  calculation 
(Equation  2) : 


XLIM(N,C,P) 


+  C 


(4) 


where 

XLIM(N,C,P)  =  nutrient  limitation  function  for  nitrogen,  carbon,  and 
phosphorus,  unitless 

C  =  concentration  of  respective  nutrient  in  the  water 

_3 

column  (N,  C,  P)  or  sediment  (N,  P) ,  g  m 
K^2  =  user-specified  half-saturation  coefficient  for  the 
respective  nutrient,  g  m  ^ 

The  limiting  nutrient  is  defined  in  this  context  as  the  one  giving  the  minimum 
value  of  Equation  4. 

11.  Many  nutrients  used  by  freshwater  submersed  macrophytes,  including 
both  nitrogen  and  phosphorus,  are  obtained  primarily  through  the  roots  from 
sediment  (Best  and  Mantai  1978;  Bole  and  Allan  1978;  Carignan  and  Kalff  1980; 
DeMarte  and  Hartman  1974;  Nichols  and  Kinney  1976).  CE-QUAL-R1  has 
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APPENDIX  A:  MACROPHYTE  MODEL  STAND-ALONE  VERSION 


Introduction 

1.  A  stand-alone  version  of  the  macrophyte  model  was  developed  to 

verify  and  validate  several  of  the  recommended  process  equations  for  a  single 
model  layer.  This  appendix  provides  a  list  of  the  state  variable  equations  * 

used  in  this  version  of  the  model.  Seven  compartments  are  represented  by  the 

model,  including  macrophytes,  dissolved  oxygen,  particulate  organic  matter 
(POM),  dissolved  organic  matter  (DOM),  phosphorus,  nitrogen,  and  organic  sedi¬ 
ment.  The  individual  process  equations  which  together  comprise  the  state 
variable  equations  are  also  described  herein.  A  parameter  list  (Table  Al) 
describes  each  of  the  parameters  used  in  the  process  equations  and  the  values 
used  in  running  the  stand-alone  version. 

2.  The  macrophyte  process  equations  correspond  to  those  given  in  the 
main  body  of  this  report  (although  several  variable  names  have  been  changed  in 
this  version  of  the  model).  Equations  for  the  other  six  state  variables  con¬ 
tain  terms  reflecting  the  impacts  of  macrophyte  processes  on  other  components 
of  aquatic  ecosystems.  This  stand-alone  version  of  the  model  is  appropriate 
for  implementation  on  a  microcomputer. 

3.  There  are  some  differences  between  the  stand-alone  version  of  the  ^ 

model  and  that  recommended  for  CE-QUAL-R1.  For  the  stand-alone  version, 

(a)  it  was  assumed  that  macrophyte  production  was  not  nutrient  limited, 

(b)  contributions  to  nutrients  from  macrophyte  respiration  were  not  included, 
and  (c)  contributions  to  nutrients  from  macrophyte  nonpredatory  mortality  are 
included.  Additionally,  CE-QUAL-R1  does  not  include  harvesting  as  the  stand¬ 
alone  version  does. 
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Recommendations 


31.  It  is  recommended  that  the  model  for  submersed  aquatic  plants 
described  in  this  report  be  incorporated  in  the  CE-QUAL-R1  model  with  due  con¬ 
sideration  of  the  following  points: 

a.  The  light  response  function  should  permit  representation  of  pho¬ 
toinhibition  (this  same  algorithm  should  be  used  for  algae  in 
CE-QUAL-R1) . 

b.  Because  nutrients  are  an  explicit  part  of  the  photosynthesis 
algorithm,  limitation  should  be  based  on  the  Monod  function  for 
the  nutrient  shown  to  be  limiting  using  threshold  ratios. 

c.  The  spatial  relationships  of  the  rooted  zone  of  macrophytes  to 
the  model  layers  should  be  accounted  for  based  on  the  intersec¬ 
tion  of  model  layers  with  the  reservoir  bottom,  creating  a  two- 
dimensional  array  of  cells  for  macrophyte  computations;  the 
macrophytes  should  be  apportioned  into  the  vertical  layers  based 
on  cell-by-cell  computations  and  a  comparison  with  a  user- 
specified  maximum  macrophyte  density  in  each  cell;  this  algo¬ 
rithm  can  also  be  used  to  determine  the  biomass  of  macrophytes 
cut  by  a  mechanical  harvester  set  at  a  particular  depth. 

d.  Chemical  control  can  be  modeled  using  dose-response 
relationships. 
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Depending  on  how  Che  chemical  control  program  is  implemented,  Che  macrophyte 
mass  killed  can  be  transferred  as  appropriate  to  other  model  compartments 
(detritus,  sediment,  dissolved  organic  matter). 

Process  Validation 


27.  Select  process  equations  included  in  the  proposed  macrophyte  sub¬ 
model  have  been  validated  based  on  experimental  results  from  the  literature 
and  published  experimental  results  performed  at  the  US  Army  Engineer  Waterways 
Experiment  Station  by  Dr.  John  Barko  and  colleagues.  Data  on  two  macrophyte 
species  of  particular  interest  to  the  Corps  were  used  in  this  validation  pro¬ 
cedure,  M.  spicatum  and  H.  verticillata .  Results  of  validating  several  spe¬ 
cific  equations  in  the  macrophyte  model  are  discussed  in  the  following 
paragraphs . 

28.  The  equation  used  to  represent  the  photosynthetic  light  response  is 
that  of  Steele  (1962)  (see  Equation  5  and  Appendix  B) .  Figures  3  and  4  demon¬ 
strate  that  this  equation  fits  experimental  data  from  Van,  Haller,  and  Bowes 
(1976)  for  M.  spicatum  and  from  Barko  et  al.  (1980)  for  H.  verticillata .  The 

parameter  PISAT,  which  describes  the  saturating  light  intensity  for  photosyn- 

-2  -1 

thesis,  was  set  at  112  and  196  kcal  m  hr  ,  respectively,  for  M.  spicatum 
and  H.  verticillata  (Appendix  B) .  Photoinhibition  at  high  light  intensities 
can  also  be  predicted  using  this  equation.  Although  this  type  of  response  of 
these  two  species  to  high  light  intensities  has  not  been  observed,  other  spe¬ 
cies  demonstrate  photoinhibition  which  could  be  significant  during  reservoir 
drawdown . 

29.  The  effect  of  temperature  on  photosynthesis  is  represented  using  the 
equation  of  Thornton  and  Lessem  (1978)  (Equation  3).  Validation  of  this  equa¬ 
tion  for  H.  verticillata ,  based  on  results  of  Barko  et  al.  (1980),  is  demon¬ 
strated  in  Figure  5.  The  parameter  values  used  in  this  equation  are  as  fol¬ 
lows:  T1  -  10°C,  T2  -  20°C,  T3  -  24°C,  T4  -  32°C,  K1  -  0.01,  K2  -  0.98, 

K3  ”  0.98,  and  K4  -  0.30  (Appendices  A  and  B) . 

30.  Validation  of  the  equation  representing  dark  respiration  (Equa¬ 
tion  6)  is  represented  in  Figure  6  for  H.  verticillata.  The  parameter  values 
used  are  as  follows:  T1  -  5°C,  T2  *  25°C,  K1  -  0.01,  and  K2  *  0.98 
(Appendix  B) . 
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The  index  J  ranges  from  1  (top  layer)  up  to  a  user-specified  value  indicating 
the  maximum  number  of  layers  in  which  macrophytes  can  occur  (actually,  the 
maximum  rooting  depth  in  metres).  If  all  the  mass  in  that  column  can  be  con¬ 
tained  in  the  bottommost  cell,  it  is  placed  there.  Otherwise,  Equation  9  is 
iterated  (i.e.,  the  value  of  J  is  increased  sequentially)  until  the  calculated 
total  macrophyte  mass  for  that  column  is  apportioned  among  cells  in  that  col¬ 
umn,  such  that  the  mass  in  each  cell  is  less  than  or  equal  to  the  maximum  cal¬ 
culated  with  Equation  9.  The  total  macrophyte  mass  is  then  calculated  for 
each  model  layer  by  summation,  and  for  the  entire  reservoir. 

Management  Control  Processes 


25.  In  addition  to  ecological  processes,  the  model  cat  also  simulate 
management  control  processes  including  mechanical  harvesting  and  chemical  con¬ 
trol  of  the  plants.  Macrophyte  mass  removed  by  mechanical  harvesting  is  a 
function  of  plant  rooting  depth  and  mass  density  as  well  as  the  cutting  depth 
of  the  mechanical  harvester.  Having  determined  macrophyte  biomass  in  each 
model  layer,  the  amount  cut  (MBIOCUT)  by  a  mechanical  harvester  set  at  a  par¬ 
ticular  cutting  depth  (CUTZ)  can  be  calculated  by  summation.  If  the  cutting 
depth  falls  between  layer  boundaries,  then  an  appropriate  fraction  of  the 
macrophyte  mass  in  that  layer  can  be  removed  since  mass  is  assumed  to  be  dis¬ 
tributed  homogeneously  within  layers. 

26.  Chemical  control  is  a  function  of  the  following  dose-response  curve 
for  the  herbicide  used: 


MCHEM(I)  -  MACRO ( I )  *  CHEM/(LC50  +  CHEM) 


where 


MCHEM(I)  ■  macrophyte  biomass  killed  in  layer  I,  g 
MACRO(I)  ■  total  macrophyte  biomass  in  layer  I,  g 

CHEM  •*  user-specified  ambient  environmental  concentration  of 
herbicide  applied,  yg  l 

LC50  ”  user-specified  herbicide  concentration  which  will  kill 
50  percent  of  the  macrophytes,  yg  i  1 
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the  Ith  layer.  The  actual  volume  of  the  Ith  layer,  DVOL(I),  is  calculated  as 
the  difference  between  V0L(I+1)  and  VOL(I).  Both  volume  and  area  are  typi¬ 
cally  represented  as  power  functions  of  elevation. 

23.  Using  this  scheme,  a  series  of  vertical  segments  or  columns  can  be 
superimposed  at  the  points  at  which  boundary  layers  intersect  the  reservoir 
bottom  (Figure  2b) ,  creating  a  series  of  two-dimensional  cells  for  macrophyte 
computations  (Figure  2c).  To  simplify  the  computational  sequence,  these  cells 
are  numbered  from  the  reservoir  surface  down,  and  from  upstream  toward  the 
dam.  A  given  cell  is  indexed  (i,j)  with  i  referring  to  row  position  and  j  to 
column.  Because  each  of  the  layers  in  the  model  representation  of  a  reservoir 
is  extremely  long  and  thin,  the  bottom  surface  area  in  which  macrophytes  root 
can  be  approximated  as  the  difference  AREA(I+1)  -  AREA(I) .  Similarly,  the 
volume  of  each  computational  cell  can  be  approximated  as  this  bottom  surface 
area  times  the  thickness  (SDZ)  of  the  layer  in  which  that  cell  occurs.  These 
bottom  areas  and  cell  volumes  are  used  in  macrophyte  computations  as  described 
in  the  following  paragraph. 

24.  Macrophytes  are  associated  with  the  bottom  sediments  in  which  they 
are  rooted  and  with  the  overlying  water  column.  In  order  to  determine  how 
macrophyte  mass  is  apportioned  among  the  cells  in  a  given  vertical  column,  the 
assumption  is  made  that  the  volumetric  density  of  macrophyte  dry  mass  cannot 

_3 

exceed  a  user-specified  maximum  value  (PLDENS,  g  m  ).  At  each  model  time 
step,  the  macrophyte  differential  equation  (Equation  1)  Is  solved  on  a  cell- 
by-cell  basis  using  a  simple  Euler  procedure  and  the  mass  Is  calculated  at  the 
previous  time  step  as  an  initial  value .  Then  macrophyte  mass  is  summed  over 
all  cells  in  a  given  column.  Beginning  with  the  bottommost  cell  (i.e.,  the 
one  nearest  the  sediment) ,  this  summed  mass  is  apportioned  among  cells  by  com¬ 
paring  it  with  the  maximum  mass  which  each  cell  can  contain.  For  cells  in  the 
Ith  column,  this  maximum  is  calculated  as 

DATA(J,I)  -  PLDENS  *  SDZ(J)  *  (AREA(I+1)  -  AREA ( I ) )  (9) 


where 

DATA(J,I) 

PLDENS 

SDZ(J) 

AREA(I) 


maximum  macrophyte  mass  which  can  be  contained  In  the 
cell  in  layer  J  and  column  1,  g 

user-specified  maximum  macrophyte  volumetric  density,  g  m 
thickness  of  Jth  model  layer,  m 
bottom  surface  area  at  layer  I,  m2 
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coefficients  are  involved.  In  a  similar  manner,  grazing  represents  a  direct 
transfer  of  mass  to  fish,  without  conversion.  As  a  consequence  of  nonpreda- 
tory  mortality,  macrophyte  biomass  is  transferred  to  dissolved  organic  matter, 
detritus,  and  sediment  compartments.  The  "dead"  biomass  is  apportioned 
between  the  three  receiver  compartments  based  on  user-specified  coefficients. 

20.  Included  in  Appendix  A  is  a  stand-alone  version  of  the  macrophyte 
model  which  was  used  in  validating  the  various  process  equations  just  dis¬ 
cussed.  In  addition  to  containing  the  equations  describing  macrophyte  physio¬ 
logical  processes  (Equations  1-8) ,  this  version  of  the  model  also  contains 
equations  for  oxygen,  particulate  organic  matter,  dissolved  organic  matter, 
phosphorus,  nitrogen,  and  sediment.  This  model  thus  illustrates  the  way  in 
which  macrophyte  terms  enter  into  equations  for  other  water  quality  constitu¬ 
ents  included  in  CE-QUAL-R1.  In  Appendix  B,  representative  values  for  the 
parameters  included  in  Equations  1-8  of  the  macrophyte  model  (as  defined  in 
Appendix  A)  are  listed,  based  on  research  on  two  macrophyte  species  of  par¬ 
ticular  interest,  Myriophyllum  spioatum  and  Hydrilla  verticillata .  CE-QUAL- 
Rl-related  parameters  and  coefficient  values  are  also  listed  in  Appendix  B. 


Spatial  Relationships 


21.  In  order  to  describe  vertical  growth  of  macrophytes  in  a  one¬ 
dimensional,  variable-layer  model  like  CE-QUAL-R1,  it  was  necessary  to  devise 
a  means  of  geometrically  segmenting  the  model  into  a  matrix  of  rows  (layers) 
and  columns.  This  matrix  defines  the  volume  of  each  segment  and  the  proximity 
of  one  segment  to  another.  A  description  of  how  the  matrix  can  be  incor¬ 
porated  into  the  CE-QUAL-R1  model  follows. 

22.  CE-QUAL-R1  is  a  one-dimensional  model  with  multiple  layers.  Thermal 
energy  and  materials  are  assumed  to  be  uniformly  distributed  within  each  model 
layer.  Reservoir  morphometry  is  represented  in  the  model  by  a  variable-layer 
approach  (i.e.,  layer  dimensions  vary  over  time  based  on  inflows  and  outflows 
and  on  user-specified  morphometric  relationships  of  area  and  volume  to  eleva¬ 
tion  above  the  reservoir  bottom).  Relationships  among  elevation,  area,  and 
volume  are  depicted  in  Figure  2a.  A  given  layer  (numbered  I,  from  the  bottom 
up)  is  specified  as  being  Z(I)  metres  above  the  bottom  and  SDZ(I)  metres 
thick.  The  area  of  the  Ith  layer,  AREA(I) ,  is  defined  at  the  lower  boundary 
of  that  layer.  A  volume,  VOL(I),  is  also  defined  up  to  the  lower  surface  of 
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calculated  as  the  product  of  the  two  temperature  limitation  functions,  RMULT1 
and  RMULT2  (Equation  3),  times  a  user-specified  maximum  fish  grazing  rate, 
times  a  Monod  function  similar  in  form  to  Equation  4.  In  this  fish-grazing 
limitation  function,  the  role  of  C  (Equation  4)  is  played  by  the  sum,  over 
all  types  of  food  (Including  macrophytes)  ingested  by  fish,  of  products  of  a 
user-specified  preference  factor  for  that  food  type  and  the  concentration  of 
that  food  type.  For  this  grazing  function,  E<Iuation  4)  would  again 

be  a  user-specified  half-saturation  coefficient  for  fish  grazing.  The  reader 
should  consult  the  CE-QUAL-R1  User's  Manual  (Environmental  Laboratory  1982) 
for  further  details.  An  additional  preference  factor  would  need  to  be  in¬ 
cluded  in  the  model,  specifying  the  fractional  preference  of  fish  for 
macrophytes. 

Interactions  with  other 
compartments  in  CE-QUAL-R1 

18.  As  depicted  in  Figure  1,  those  macrophyte  processes  discussed  above 
also  impact  a  variety  of  other  compartments  in  CE-QUAL-Rl.  Thus,  correspond¬ 
ing  to  the  process  equations  given  above  (Equations  1-8) ,  terms  will  need  to 
be  added  to  or  subtracted  from  other  equations  in  the  model.  These  terms  rep¬ 
resent  the  addition  or  removal  of  mass  to  or  from  other  compartments  in  the 
modeled  reservoir.  These  terms  will  be  briefly  described  here.  Although  the 
actual  equations  will  not  be  provided,  they  correspond  exactly  to  the  form  of 
the  equations  listed  previously. 

19.  As  a  result  of  macrophyte  photosynthetic  processes,  oxygen  is 
evolved.  This  is  modeled  as  an  "equivalent  oxygen  concentration,"  calculated 
as  the  product  at  the  gross  production  rate  of  concentration  and  a  user- 
specified  oxygen-to-biomass  stoichiometric  coefficient,  which  is  added 
directly  to  the  oxygen  differential  equation.  Similarly,  dark  respiration 
removes  oxygen.  This  removal,  a  subtraction  from  the  oxygen  equation,  is  cal¬ 
culated  as  the  product  of  the  dark  respiration  rate  of  concentration  and 
another  user-specified  stoichiometric  coefficient.  Gross  production  and 
respiration  also  result  in  the  uptake  and  release,  respectively,  of  nutrients 
(N,  P,  C)  from  and  to  the  water  column  and  sediments  (Figure  1).  These  trans¬ 
fers  are  calculated  as  the  product  of  the  production  and  respiration  rates  of 
concentration  and  user-specified  nutrient-to-biomass  stoichiometric  coeffi¬ 
cients.  Photorespiration  represents  a  direct  addition  of  mass  to  the  ammonia- 
nitrogen,  phosphorus,  and  dissolved  organic  matter  compartments;  no  conversion 
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mathematically  using  only  the  rising  limb  of  the  temperature  equation  of 
Thornton  and  Lessem  (1978)  (Equation  3): 

MRESP  =  MKhESP  *  RMULTl(T)  (6) 

where 

MRESP  =  dark  respiration  rate,  day  * 

MKRESP  =  user-specified  maximum  dark  respiration  rate,  day  * 

15.  Photorespiration.  Photorespiration  or  excretion  is  important  be¬ 
cause  it  results  in  the  phenomenon  known  as  "nutrient  pumping,"  whereby  nu¬ 
trients  are  transferred  from  bottom  sediments  to  water.  This  process  also 
increases  the  amount  of  organic  matter  dissolved  in  the  water  column.  Excre¬ 
tion  is  a  function  of  light  intensity.  Under  conditions  of  very  high  or  very 
low  light  intensities,  the  rate  of  extracellular  release  increases.  Mathe¬ 
matically  this  is  represented  as 

MEXCR  =  (1  -  XLIML)  *  MKEXCR  (7) 

where 

MEXCR  =  excretion  rate,  day 
MKEXCR  =  user-specified  maximum  excretion  rate,  day  1 

16.  Nonpredatory  mortality.  Nonpredatory  mortality  is  temperature- 
dependent  when  the  change  in  temperature  (increase  or  decrease)  over  a  7-day 
period  exceeds  a  critical  maximum  temperature  TMPMAX  .  Therefore,  if 
|TMPTUR(1)  -  TMPTUR ( 7) |  >  TMPMAX: 

MMORT  -  MKMORT  (8) 


where 

TMPTUR ( 1 )  and  TMPTUR(7)  ■  water  temperature  over  7-day  period,  °C 
TMPMAX  -  maximum  temperature  change,  ®C 
MMORT  -  nonpredatory  mortality  rate,  day  * 

MKMORT  ■  user-specified  maximum  nonpredatory  mortality 
rate,  day  * 

17.  Grazing.  Grazing  of  macrophytes  by  fish  is  modeled  with  the  same 
type  of  grazing  function  as  used  in  CE-QUAL-R1.  Thus,  the  grazing  rate  is 
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compartments  representing  sediment  nitrogen  and  phosphorus;  therefore,  limita¬ 
tion  of  nutrients  obtained  through  the  roots  can  occur,  although  this  is  rare 
in  nature.  This  process  is  most  important  in  allowing  "nutrient  pumping"  from 
the  sediments  into  the  water  column. 

12.  In  some  cases  where  nutrient  concentrations  in  the  water  are  high, 
it  becomes  advantageous  for  the  plant  to  draw  nutrients  from  the  water  column. 
In  water  with  a  phosphorus  concentration  of  2.0  mg  l  ,  characteristic  of 
eutrophic  reservoirs,  Myriophyllum  spieatum  took  phosphorus  from  the  water 
column  (Bole  and  Allan  1978) .  This  is  modeled  using  a  species-specific 
parameter  to  indicate  the  water  concentration  above  which  nutrients  are  taken 
from  the  water  column.  Whenever  the  water  column  concentration  of  nitrogen  or 
phosphorus  equals  or  exceeds  this  user-specified  concentration,  it  is  the 
water  concentration  of  that  nutrient  which  is  entered  into  the  Monod  equation 
(Equation  4).  Otherwise,  it  is  the  sediment  concentration  of  nitrogen  or 
phosphorus  which  is  used  in  Equation  4. 

13.  Light  limitation  is  represented  using  Steele's  equation  (1962): 

/0.5  *  SWSA)  /0.5  *  SWSA \"1  ... 

XLIML  -  plsAT  J  exp  |^1  -  y  plsAT  jj  (5) 


where 

-2  -1 

SWSA  -  average  irradiance  for  a  specific  model  layer,  kcal  m  hr 
(calculated  in  Subroutine  HEAT  in  CE-QUAL-R1) 

PISAT  -  user-specified  irradiance  level  at  which  the  photosynthetic  rate 

-2  -1 

is  saturated  (i.e.,  occurs  at  maximum  rate),  kcal  m  hr 

The  coefficient  value  0.5  is  used  in  Steele's  equation  to  represent  the  frac¬ 
tion  of  total  irradiance  that  is  photosynthetically  active  radiation  (PAR). 

PAR  is  in  the  range  of  400  to  700  nm.  Steele's  equation  can  predict  photoin¬ 
hibition  of  photosynthesis  at  high  light  intensities,  above  the  level  speci¬ 
fied  by  PISAT.  Solar  radiation  is  distributed  vertically  in  the  water  column 
in  CE-QUAL-R1  based  upon  the  extinction  coefficient  for  water.  Light  is  also 
attenuated  by  self-shading  by  algae,  zooplankton,  detritus,  and  suspended 
solids.  An  additional  self-shading  coefficient  should  be  included  in  the 
model  to  account  for  the  effect  of  macrophyte  biomass  on  light  attenuation. 

14.  Dark  respiration.  Dark  respiration  is  a  function  of  temperature. 

As  with  other  respiratory  rates  in  CE-QUAL-RI,  it  is  represented 
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State  Variable  Equations* 


Macrophyte 

MACRO  =  MPROD  -  MRESP  -  MEXCR  -  MMORT  -  MHVST 

macrophyte  =  photosynthesis  -  dark  respiration  -  excretion  (photo¬ 
respiration)  -  mortality  -  harvesting 

Oxygen 

OXYGEN  =  OTST  +  OMAC  -  ANIT  -  OPDK  -  ODDK  -  OSDR 

Oxygen  =  oxygen  saturation  +  contribution  from  macrophytes 

-  equivalent  loss  from  nitrogen  decay  -  equivalent  loss  from  POM  decay 

-  equivalent  loss  from  DOM  decay  -  equivalent  loss  from  sediment  decay 
Particulate  organic  matter 

POM  =  PMAC  -  PDK  -  PSTL 

POM  =  contribution  from  macrophyte  mortality  and  harvesting 

-  loss  from  POM  decay  -  loss  from  settling 
Dissolved  organic  matter 

DOM  =  DMAC  +  DEXCR  +  DDK  -  DBAC 

DOM  =  contribution  from  macrophyte  mortality  and  harvesting  +  contribu¬ 
tion  from  macrophyte  excretion  +  contribution  from  POM  decay  -  loss  from 
bacterial  respiration 
Phosphorus  (water  column) 

P04  »  FMAC  +  FDK  +  FEXCR  -  FSINK 

P04  =  contribution  from  macrophyte  mortality  and  harvesting  +  contribu¬ 
tion  from  decay  of  POM  and  sediments  +  contribution  from  macrophyte 
excretion  -  loss  to  algal  production 


*  Each  equation  represents  the  Jime  rate  of  change  of  the  state  variable  for 
a  model  layer.  The  units  of  MACRO  are  grams  per  square  metre  per  day  per 
model  layer.  The  units  of  all  other  state  variables  are  grams  per  square 
metre  per  day  per  metre  of  model  layer. 
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Nitrogen  (water  column) 

N  -  NMAC  +  NDK  +  NEXCR  -  NS INK 

N  «  contribution  from  macrophyte  mortality  and  harvesting 
+  contribution  from  decay  of  POM  and  sediments 

+  contribution  from  macrophyte  excretion  -  loss  to  algal  production 
Organic  sediment 

SED  -  SMAC  -  SDK 

SED  =  contribution  from  macrophyte  mortality  and  harvesting 
-  loss  from  sediment  decay 

Process  Equations 


Macrophyte 


MPROD  -  PMAX  *  RMULTl(T)  *  RMULT2 (T)  *  LIGHT  *  MACRO 


where 

PMAX  «  maximum  photosynthetic  rate,  day  1 

RMULTl(T)  »  temperature  limitation  function,  unitless 

RMULT2(T)  ■  temperature  limitation  function,  unitless 

T  ■  ambient  water  temperature,  °C 

LIGHT  ■  light  limitation  function,  unitless 

-2 

MACRO  •  macrophyte  biomass,  g  m 


LIGHT 


where 

e  -  extinction  coefficient 

Z2  -  depth  at  the  bottom  of  the  simulated  section,  m 
Zi  “  depth  at  the  top  of  the  simulated  section,  m 
10  >  irradiance  at  the  water  surface,  kcal  m  sec 

-2 

ISAT  -  saturating  irradiance  for  photosynthesis,  kcal  m  sec 


MRESP  -  KRESP  *  RMULT1 (T)  *  MACRO 


where 

KRESP  =  user-specified  maximum  respiration  rate,  g  g  '  day 
MEXCR  “  KEXCR  *  (1-LIGHT)  *  MACRO 

where 

KEXCR  -  user-specified  maximum  excretion  rate,  g  g  1  day 
If,  |TMPTUR(1)  -  TMPTUR(7) j  is  greater  than  TMPMAX,  then 

MMORT  -  KMORT  *  MACRO 


where 

KMORT  »  nonpredatory  mortality  rate,  g  g  1  day  * 

TMPMAX  •>  critical  maximum  temperature  difference  over  a 
7-day  period,  °C 

TMPTUR ( 1 )  and  TMPTUR(7)  -  water  temperatures  over  a  7-day  period,  °C 

MHVST  -  CHEM  *  MACRO 


where 

CHEM  =  rate  of  die-off  of  macrophyte  dependent  upon  type  of  chemical 
used,  g  g  '  day  1 

NOTE:  Mechanical  harvesting  is  calculated  outside  the  differential 
equation  as  follows: 


MWH  +  MACRO  -  MHT 
Z  -  MHT  -  TPLT 
CUTZ  -  TPLT  -  MCUT 
MWH  i  MCUT  -  MBIOCUT 


where 


MWH  -  species-specific  weight-to-height  ratio,  g  n‘ 
MHT  -  macrophyte  height,  m 
2  «  depth  of  water  column,  m 
TPLT  »  top  of  plant,  m 

CUTZ  ••  cutting  depth  of  mechanical  harvester,  m 
MCUT  -  amount  of  macrophyte  cut,  m 
MBIOCUT  *  biomass  of  macrophyte  cut,  g  m-4 


Oxygen 


OTST  =  (14.6  *  exp(-(0. 027767  -  0.00027  *  T 
+  0.000002  *  T  *  T)  *  T))  *  Z 
OMAC  =  (0MACEQ1  *  MPR0D)  -  (0MACEQ2  *  MRESP) 
AN IT  -  ONEQ  *  NMAC 
OPDK  «*  OPEQ  *  PDK 
ODDK  =  ODEQ  *  DDK 
OSDK  »  OSEQ  *  SDK 


Particulate  organic  matter 

PMAC  -  (MMORT  *  Ml)  +  (MHVST  *  HI) 

PDK  -  RPOM  *  POM  *  RMULTl(T) 

PSTL  «  (PMSTL  *  MMORT)  +  (PHSTL  *  MHVST) 
KPOM  -  0.01192  *  1/NTC(2)  +  0.00672 


Dissolved  organic  matter 

DMAC  =  (MMORT  *  M2)  +  (MHVST  *  H2) 
DEXCR  -  MEXCR  *  E2 
DDK  -  PDK  *  P2 

DBAC  -  KD0M  *  DOM  *  RMULTl(T) 

KDOM  -  0.024  *  1/NTC(3)  +  0.0192 

Phosphorus 

FMAC  -  (MMORT  *  M3)  +  (MHVST  *  H3) 
FDK  -  (PDK  *  P3)  +  (SDK  *  S3) 
FEXCR  -  MEXCR  *  E3 
FSINK  “  photoplankton  biomass  *  FRS 
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Nitrogen 


NMAC  -  (MMORT  *  M4)  +  (MHVST  *  H4) 
NDK  -  (PDK  *  P4)  +  (SDK  *  S4) 
NEXCR  -  MEXCR  *  E4 
NS INK  ■  photoplankton  biomass  *  NRS 


Sediments 

SMAC  -  (MMORT  *  M5)  +  (MHVST  *  H5) 
SDK  -  KSED  *  SED  *  RMULTl(T) 

KSED  -  0.00519  *  1/NTC(4)  ♦  0.00346 


Stand-Alone  Varalon  Macrophyte  Model  Parameter  List 


Table  A1  (Continued) 


photosynthesis  and  0.98  Calibrated 


Tabic  A1  (Concluded) 


APPENDIX  B:  MACROPHYTE  MODEL  PARAMETER  LIST 


Tabulated  In  Table  BI  In  this  Appendix  are  values  for  specific  parameters 
included  in  the  state  variable  and  process  equations  vhlch  comprise  the  macro¬ 
phyte  model  proposed  in  the  main  body  of  this  report  (as  intended  for  inclu¬ 
sion  in  CE-QUAL-R1) .  These  values  were  either  derived  from  published  litera¬ 
ture  sources  or  established  in  the  process  validation  studies  described 
earlier.  Most  values  tabulated  here  apply  to  one  or  two  macrophyte  species  of 
interest,  Myriophyllum  spicatwn  or  Hydrilla  verticil  lata. 
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Parameter  Description  _ Species _  _ Value _  Converted  Value _ Reference 
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